Transport Statistics of Bistable Systems 
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We consider the transport statistics of classical bistable systems driven by noise. The stochastic 
path integral formalism is used to investigate the dynamics and distribution of transmitted charge. 
Switching rates between the two stable states are found from an instanton calculation, leading to 
an effective two-state system on a long time scale. In the bistable current range, the telegraph noise 
dominates the distribution, whose logarithm is found to be universally described by a tilted ellipse. 
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Bistable systems are ubiquitous in nature. They occur 
in solid state [l|, |3| , magnetic |3 , and biological systems 
|j| . The behavior of bistable systems may be character- 
ized by random switching between two stable states, in- 
duced by either thermal fluctuations, or non-equilibrium 
noise. The presence of bistability is usually deduced by 
measuring the current passing through the system. The 
switching of the internal state of the system alters the 
resistance of the sample, causing the measured current 
to randomly switch between two distinct average values. 
The noise of this switching current is called random tele- 
graph noise |5J, a process studied by Machlup already in 
the 1950s. An important feature of this process is that 
its zero frequency noise power diverges as the switching 
rates tend to zero. 

Fluctuations in bistable electronic devices, such as tun- 
nel diodes, were theoretically investigated by Landauer 
[g, and later by Hanggi and Thomas |7|. Since then, 
bistable fluctuations have been observed in microstruc- 
tures by different experimental groups [S|. The switch- 
ing times were found to be well fit by an exponential of 
the applied bias. More recently, similar physics was dis- 
covered in mesoscopic resonant tunneling wells near an 
instability J2|. In this Letter, we examine the properties 
of bistable systems from the novel point of view of trans- 
port statistics. This approach gives not only the average 
current, and zero- frequency noise [lOJ, but all higher cur- 
rent cumulants as well (also called full counting statistics 
in mesoscopic physics |ll|). For many systems, high or- 
der cumulants may be computed via the cascade method 
13, which uses the idea that the noise itself has noise 
13. I13L Il4| . and that feedback effects of the lower or- 
der cumulants may be taken into account in a perturba- 
tive manner [isL llq. In contrast, for unstable systems, 
non-linearities enter the transport in an essentially non- 
perturbative fashion. Here the situation is reversed: All 
higher cumulants of the noise sources must be accounted 
for to correctly find the switching rates which govern even 
the average current. 

In order to successfully describe the physics of such an 
instability, we apply the Stochastic Path Integral (SPI), 
a non-perturbative method that was recently introduced 
llM llol to evaluate the statistical fluctuations in a clas- 




FIG. 1: A schematic of the average current flowing through 
the left and right connector as a function of the charge in 
the central node. While linear connectors (the dashed line) 
have only one crossing point, nonlinear devices may have three 
crossing points, two stable, Qi,2, and one unstable, Qu- Inset: 
The stochastic network. The direction of current flow through 
the network is indicated. 



sical stochastic network. The network is composed of 
nodes containing effectively continuous charge represent- 
ing any conserved quantity. The charge is transferred 
through connectors whose isolated transport properties 
are given. In order to emphasize the universality of the 
results and focus on the essential physics, we use an ab- 
stract description which can represent many of the phys- 
ical systems mentioned thus far, such as the tunneling 
diode, or more general S or Z-type instabilities ,2j . 

We now briefly describe our main results. We let the 
average left current have a common S-type nonlinear de- 
pendence on the charge on the nodes (see Fig. QJ , and 
demonstrate bistability. The connector noise induces 
switching between the stable states with the rates ri^2, 
found via an instanton calculation, Eq. Q, to be ex- 
ponentially small with the instability charge scale. On 
the time scale of order rj~2 or longer, the system can 
be viewed as being an effectively two-state system. Af- 
ter this point, results H12I15|I hold independently of how 
the rates ri^2 are obtained. The requirement of contin- 
uous charge may be relaxed, and the results apply to 
any telegraph process [ll |2) S 14 with the rates as input 
parameters. The distribution of the transmitted charge 



is dominated by one of the stable states in the extreme 
small or large current ranges. In between, the distribu- 
tion has a comparatively flat region dominated by the 
telegraph process, Eq. l(T^ . Given ri_2, it is universal: 
The logarithm of the distribution is a tilted ellipse. 

The model.-The bistable network we consider con- 
sists of one "conserving" node, C, storing the general- 
ized charge Q. A typical non-equilibrium fluctuation de- 
cays on a time scale tq. This node is connected to two 
"absorbing" nodes, L and R, which inject and absorb 
currents II and Ir through the connectors between the 
nodes (see inset of Fig. QJ. The sources of noise are the 
currents Il,Ri which are Markovian random variables (af- 
ter a time scale tq <C tc) with known statistics given by 
the current cumulant generating functions, 
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where the cumulants ((la)) arc functions of the charge in 
the central node, Q [llf. Charge conservation, Q = Il~ 
Ir, leads to the slow dynamics of Q which couples back 
into the dynamics of Il,r by changing the conditions 
for transport through the connectors, thus altering the 
statistics of the net current / = (Jl + ^fi)/2. The cumu- 
lants ((Q")) of the transmitted charge Q{t) — Jp dt'I{t') 
can be obtained with diagrammatic perturbation theory 
|l2Lll6l | in the stable regime by using the SPI formalism. 
However, the power of the SPI method is that it also al- 
lows for a non-perturbative analysis of unstable systems. 
The SPI represents the time evolution of the system 
while monitoring all fluctuations 15, 16] : 



u[xA 



VQVX exp 



dt'[-\Q + H{Q,\x)]], 



with fixed initial charge Q{Q) and the final condition 
\{t) = 0. The generating function S'(x, t) for the cumu- 
lants of the transmitted charge ((Q")) = d"5(x,0lx=o 
is given by S' = ln[/[x,i]. The path integral has a 
canonical form, with the "Hamiltonian" H given by 

H(Q, A, x) = Hl{\ + x/2) + Hr{-\ + x/2), (3) 

where Hl.r are the current cumulant generators (^. The 
counting variable A is a Lagrange multiplier conserving 
the charge on the node C . A large parameter, the num- 
ber of elementary charges participating in the transport, 
allows the saddle-point evaluation of the SPI |5J leading 
to the canonical equations of motion 



Q = dH/d\ A = -dH/dQ. 



(4) 



The solution of these equations should be substituted 

into Eq. (0) in order to obtain the generating function S. 

The bistability.-We first assume for simplicity that the 

sources /„, a = L,R are Gaussian with average current 




FIG. 2: The stability profile. The Hamiltonian is given as a 
phase space density plot. The system has three fixed points, 
two stable ai Q — Qi,Q2 and one unstable at Q = Q„. For 
long time, the system is forced to the zero energy lines in 
black, Ao,Ain. A decay from the unstable state may happen 
with zero action, but to hop from either stable state to the 
unstable state requires finite action, the area underneath the 
curves shown. 

Ia{Q) = (la) and noise power Fa{Q) — (({laf)), and 
then generalize to arbitrary noise generators (^. Fixing 
X = 0, we find that the unconditional dynamics of the 
system (where the transmitted charge is unmonitored) is 
governed by the Hamiltonian 



Ho - {II - Ib.)X + (1/2) (Fi + Fr)\'' 



(5) 



Therefore, the stationary state Q = -^ = is [according 
to Eq. (gl)] given by A = and Il{Q) = Ir{Q), the con- 
servation of average current. In the ohmic regime, the 
average current is a linear function of the charge in the 
node (see Fig. ^| , so the currents intersect only at one 
point. However, if we allow a nonlinear I-V, topological 
distortions of the curves allow three crossing points, giv- 
ing three stationary solutions. We will now show that 
while the outer two intersections are stable to perturba- 
tions, the central intersection is unstable. 

First, we note that because the Hamiltonian (|SJ) does 
not explicitly depend on time, it is an integral of mo- 
tion: Ho{Q,X) = E, where the "energy" determines a 
trajectory in phase space. On a time scale longer than 
the relaxation time, t ^ tc, the dissipative dynamics of 
the system projects a Q-distributed initial state onto the 
zero energy lines, Ho{Q,X) — 0. According to Eq. ^, 
there are at least two such lines: The Aq line with A = 0, 
and the "instanton" line Ain, given by 



Ai„(0) = -2(/l - lR.)/iFL + Fr). 



(6) 



These lines, shown in Fig.|21 connect the three stationary 
points, Qi, Qu, and Q2- From Eqs. Q, for the line Aq 



we have Q = II — Ir, while for the Am hne the sign is 
opposite, Q = Ir — II- Then the simple investigation of 
Fig. n shows the propagation direction. 

According to Eq. Q, the action on the zero energy 
lines is given by 5* = — J X dQ. On the Aq line, the action 
vanishes, which makes the dynamics along this line most 
probable. This dynamics describes the relaxation of the 
initial state to the points Qi and Q2 [see Fig. (01], mak- 
ing them the stable stationary points of the system, and 
Qu the unstable point. The propagation along the Ain 
line away from stable points generates a non-zero action, 
and therefore it has an exponentially small probability. 
Occasionally, the system switches between points Qi and 
Q2 on the lines Aq, Ain with the action S = —Ai^2, 



Qu 

Al,2= / dgAin(Q), 

'Q1.2 



(7) 



given in Fig. [3 by the area between the hnes Aq and Ain. 

For an estimate, we consider the simplest (symmetric 
cubic) current difference, II — Ir = —aQ[Q^ — (AQ)^], 
and constant noise, Fl + Fjj = 2F. For the action (0 
with the instanton @, wc obtain Ai^2 ~ (a/4F)(A(5)^. 
By setting A1.2 ~ 1 and using Eqs. I@J| we estimate 
the switching time Tin ^ (AQ)'^ /F, a new time scale in 
the problem. The assumption of effectively continuous 
charge demands that AQ ^ 1, leading to Tin ^ To, so 
that the separation of time scales requirement is satisfied. 
This also implies the instability is well developed. 

The result ^ derived above for Gaussian noise with 
the instanton (jBJ holds for general noise with the Ain line 
defined by the equation 



Ho{Q,X)=Hl{X)+Hr{-X) = 0. 



(8) 



From now on we consider general noise H^ji ^18.] . 

The effective two-state Hamiltonian.-We now proceed 
with the instanton calculation [\^. We consider a time 
interval Ai much longer than the switching time Tin- On 
this time scale, the system always relaxes to one of the 
stable points (5i,2, so that it can be considered as being 
effectively a two state system where the evolution opera- 
tor (0) becomes a 2 x 2 matrix. Starting from one of the 
stable points, the system will switch to the other stable 
point with the small transition probability ri_2Ai that is 
proportional to the time interval, because the attempts 
are uncorrelated. The switching rates are given by the 
instanton (Ain-line) contribution to the path integral. 



ri.2 = wi,2exp(-Ai^2), 



(9) 



where the action Ai^2 is given by Eq. Q [23 • We note 
if the noise goes to zero, then Fi 2 vanish exponentially. 
The prefactor wi.2 is the attempt rate, which is subdom- 
inant [2l| . Finally, using the smallness of the switching 
rates, we can write a differential master equation for the 



evolution operator, U ~ HqU, where the effective Hamil- 
tonian matrix has components ±ri^2- The evolution op- 
erator describes the relaxation with the rate Ts = Fi +r2 
of any given initial state to the stationary state with con- 
stant occupation probabilities. 



Pi = r2/Fs, P2 = Fi/Fs. 



(10) 



Having described the bare (unconditional) dynamics of 
the system, we now concentrate on the transport statis- 
tics. A X contribution to the Hamiltonian ((2J) deforms the 
stability profile in Fig. |21 The stable stationary points 
acquire a non-zero energy, and thus contribute to the di- 
agonal values in the full effective Hamiltonian H. The 
rates ri^2 will also depend on x- However, the bistabil- 
ity is only important for x oc Fi^2, and therefore the x 
correction to ri^2 is exponentially small and may be ne- 
glected. Using Eq. (PJ to evaluate the contribution of the 
stable points to the action, we arrive at the result 



H 



^iii(x)-Fi r2 ^ 
, Fl H2{x)~r2, 



(11) 



The Hamiltonians i?i,2(x) are the noise generators at the 
stable points Q1.2 with the net currents Ii^2 = II = Ir. 
and noise power Fi,2 = {FlG\ + FrGI)/'{Gl + Gr)\ 
where Gl.r are the differential conductances taken at 
Qi.2- High cumulants may be found perturbatively 16]. 
The transport statistics- The Hamiltonian pi(l deter- 
mines the conditional dynamics of the effective two-state 
system on the time scale t ^ Tin [22] • It describes the 
relaxation of the system to the new conditional station- 
ary state. In an experiment on the bistable system, if 
the averaging time of the measurement is much shorter 
then the relaxation time ts = (Fi + F2)^^, the output of 
the measurement device will only reveal the noise prop- 
erties of one of the states. On the time scale i ^ rg 
the system switches many times and loses its correla- 
tion with the initial state. As a result, the transport 
becomes again Markovian, S{x) = tll{x)- The generat- 
ing function H{x) is given by the largest eigenvalue of 
the Hamiltonian (|ll|l . 

H{x) = (1/2) Y. (Hn - r„) 



(ff2-ifi-AF)V4 + FiF2, (12) 

where AF = F2 — Fi. We immediately see that the dis- 
tribution is properly normalized because 11(0) — 0. The 
average current is given by (/) — H'{0) — J2n=i 2 InPn, 
and the noise power ((/^)) = H"{0) by 



((/2)) = J2 FnPn + 2(A/)2FiF2/F|, 



(13) 



n=l,2 



where A/ = I2 — Ii, and the stationary probabilities P„ 
are given in Eq. (|10|l . The first term is the weighted 




FIG. 3: The log-distribution of transmitted charge, divided 
by ty/TiT2- Within the bistable range, the distribution is well 
fit by a tilted ellipse, connecting to the distribution dominated 
by state 1 or 2. The average current is denoted by (J). Inset: 
The same distribution, \nP[It]/t, showing the comparatively 
flat bistable distribution with the left and right asymptotics. 
ra/Ti = 4, Fi,2 = (1, 1.5) X 10-"{(/')). 



noise of the stationary states, and the second term is 
the well-known result for zero-frequency telegraph noise 
|5| which dominates the first term because the rates are 
small. The dominant contribution to the third cumulant 
is {{P)) = 6(A/)3rir2Ar/r|. The slow switching rates 
cause all cumulants to be parametrically large, making 
bistable systems a promising candidate for an experimen- 
tal measurement of the transport statistics. 

We now Fourier transform U(ix) using the stationary 
phase approximation to obtain the probability to trans- 
mit charge Q = It through the system in time t. Outside 
the bistability region, I < Ii ot I > I2 the generator in 
Eq. (fT^ may be replaced with iJi — Fi or H2 — T2, so that 
the extreme value statistics is dominated by one of the 
stable states. Within the bistable region, Ii < I < I2, 
the probability Vnil) of occupying state n — 1,2 under 
the condition that a charge It is transmitted is given by 



Vi.2{i) = ^2,1/(^1 + ^2), Gn = Vr„|/-/„|, (14) 

which generalizes (|10|l . The probability distribution of 
transmitted charge is given by 



inP{it)/t^~igi-g2f/Ai 



(15) 



and is determined by the telegraph process. The loga- 
rithm of the exact distribution (divided by the rates) is 
shown in Fig. 3 inside the bistable region and is well fit 
by a tilted ellipse, while the inset compares left and right 
asymptotics to the full distribution. 

In conclusion, we have discussed a general model of 
a noise driven bistable system. We have calculated the 
switching rates ® , transport statistics l(T^ and ((T5|l . and 
conditional occupation probabilities H14|l . The higher cu- 
mulants in Eq. Q) determine the instanton lines in Fig. 
2 [see Eq. ©], and thus make an exponentially large con- 
tribution to the rates Fi_2 tOill^- Away from the bistable 
region, the probability tails decay with a non-Gaussian 



distribution determined by the isolated statistics -ffi_2 of 
the stable states. Within the bistable current range, the 
distribution is always an ellipse whose dimensions are 
fixed by the rates alone, showing a robust universality. 

This work was supported by the SNF, and by INTAS 
(project 0014, open call 2001). 
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